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Abstract 

Recently, we pointed out that chiral transformation properties of strong penguin operators 
change in the transition from unquenched to (partially) quenched QCD. As a consequence, 
new penguin-like operators appear in the (partially) quenched theory, along with new low- 
energy constants, which should be interpreted as a quenching artifact. Here, we extend 
the analysis to the contribution of the new low-energy constants to the K° — > 7r + 7r~ 
amplitude, at leading order in chiral perturbation theory, and for arbitrary (momentum 
non-conserving) kinematics. Using these results, we provide a detailed discussion of the 
intrinsic systematic error due to this (partial) quenching artifact. We also give a simple 
recipe for the determination of the leading-order low-energy constant parameterizing the 
new operators in the case of strong LR penguins. 
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1 Introduction 



A reliable calculation of long-distance contributions to non-leptonic kaon-decay rates, and, 
in particular, to the CP- violating part parametrized by the quantity e' / e has been a long- 
standing challenge. Ideally, one would expect such calculations to be in the domain of lat- 
tice QCD, but in practice many theoretical and numerical difficulties have made progress 
in this direction rather slow. Recently however, two lattice collaborations have reported 
on numerical results for both the real and imaginary parts of AI = 1/2 and AI = 3/2 
K — > mi matrix elements with a rather satisfactory control over statistical errors [1, 2]. 
These lattice computations were done with the effective weak AS = 1 hamiltonian with 
three flavors, i.e. with the charm integrated out, and they were possible because of the use 
of lattice fermions with good chiral symmetry. Both groups reported values of e' / e which 
are non-zero, and thus consistent with the existence of direct CP violation, but of opposite 
sign (and comparable size) to the experimentally measured value. 

While statistical errors for these lattice computations seem to be reasonably under control, 
this is not the case for a large class of systematic errors, which will need to be studied further 
in the future. One source of systematic error is the use of the quenched approximation. 
In a previous paper [3] we pointed out that, in addition to the fact that quenched QCD is 
just not the same theory as full QCD, an ambiguity arises in the definition of the quenched 
version of penguin operators appearing in the AS = 1 effective weak hamiltonian. The 
ambiguity originates in the difference of the chiral transformation properties of penguin 
operators within the quenched and unquenched theories. This implies that not only will 
quenched lattice results be hampered by the fact that we do not really know whether 
quenched values of given matrix elements are close to their real- world values, but, in the 
case of penguins, they also depend on which definition of the operators is chosen, since 
more than one definition is possible. 

In fact, both lattice computations [1, 2] did not directly compute K — > tttt matrix elements, 
but K — > 7r (with M n = Mk) and K — > transition amplitudes, and used chiral pertur- 
bation theory (ChPT) to convert them into the desired K — > tttt matrix elements [4]. In 
ref. [3] we explained how the chiral properties of penguin operators change in the transition 
to the (partially) quenched theory, and how, in principle, more than one definition of a 
quenched penguin operator is possible. Using ChPT, we traced how this affects K — > ix 
and K — > matrix elements. We restricted ourselves to LR penguins {i.e. Q 5 and Qe), 
because the effects in this case already appear at leading order in ChPT, while they are 
a next-to-leading order effect for LL penguins. Because of the fact that the ambiguity is 
already present at leading order for matrix elements of Q$fi, this may be an important issue 
for e'/e (while it is expected to be less important for the real parts of K — > nn amplitudes, 
and thus the AI =1/2 rule) . 

In this paper we extend our ChPT calculations to the effect of the ambiguity on K — > Tin 
matrix elements, again to leading chiral order. This is important for two reasons. First, 
lattice computations may be done directly for K — > mr matrix elements, and their chiral 
behavior needs to be known in order to fit lattice results as a function of quark masses. 
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Since lattice computations are typically done with unphysical (i.e. energy /momentum non- 
conserving) kinematics, we present our results for the most general kinematics possible, 
both in the quenched and partially quenched cases. Second, once the complete (leading- 
order) ChPT expressions for K — > 0, K — > i\ and K — > tttt matrix elements are available, 
it is possible to give a more detailed discussion of the systematic error introduced by the 
ambiguity in the definition of quenched penguin operators. 

The paper is organized as follows. In section 2, we review the main observation of ref. [3]. 
We show how a LR penguin, which transforms in an irreducible representation (irrep) of 
SU (3)l x SU (3)r, splits into two operators in the (partially) quenched theory, with each of 
them transforming in a different irrep of the (partially) quenched chiral symmetry group. 
One of these irreps corresponds "naturally" to the single irrep of the unquenched theory, 
while the other irrep can be considered as "new," and an artifact of quenching. We give 
the ChPT realization of all relevant operators at leading and next-to-leading chiral order, 
introducing new low-energy constants (LECs) which appear in correspondence to the new 
irrep. In section 3, we present our results for K° — > 7r + 7r~ penguin matrix elements to 
leading order in ChPT, with general kinematics, and specialize these results to physical 
(i.e. energy-momentum conserving) kinematics. In section 4, we discuss different strategies 
available for using quenched lattice results to estimate real-world K — > im matrix elements, 
and give some numerical examples. We provide a simple prescription for determining the 
leading-order LEC representing the new irrep in section 5, and section 6 contains our 
conclusions. Some of this work has already been presented in ref. [5]. 



2 Review of LR penguins in (partially) quenched QCD and 
ChPT 

A lagrangian definition for partially quenched QCD can be constructed as explained in 
ref. [6] (see also ref. [7] for an alternative realization using the replica method). In addition 
to the valence quarks q vi , i — u,d, s, with masses m vi , one introduces a separate set of sea 
quarks q si , i = 1, . . . , N, with masses m si , and a set of "ghost" quarks q gi , i — u,d, s, with 
masses equal to those of the valence quarks m g i = m V i [8]. Ghost quarks are given bosonic 
statistics, such that the ghost-quark determinant cancels the valence-quark determinant, 
thus leaving only the sea-quark determinant present in the path integral. Therefore, only 
sea quarks propagate in internal loops. 

Since partially quenched QCD thus contains more flavors than unquenched QCD, its flavor 
symmetry group is larger than the QCD one. The full chiral symmetry group relevant for 
light meson physics is the graded extension of the ordinary chiral group SU(3 + N\3)l x 
SU(3 + N\3)r [6]. It is graded because part of its elements transform fermions into bosons 
and vice versa. The quenched theory, which has no sea quarks at all, corresponds to the 
special case N = [9] . 
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We consider LR penguin operators of the form 



Qpenguin = (sd) L (uU + dd + Ss)r , (2.1) 

where 

(QiQ2)l,r = QiI^Pl^ , (2.2) 
Pl,r = i^-Tls) , 

and color contractions are not specified, so that Q pen guin can represent both Q 5 and Q 6 . 
As already pointed out in ref. [3], the u, d and s fields in eq. (2.1) represent valence quarks 
in the partially quenched theory, and the penguin operator can be decomposed as 1 

3 — — — — 

Qpenguin = Stl (A^^7 M P i) str {^1^ P r) + 8t* (A^T^) Stl (Ai^^Pr) , 

Q PQS , Q PQA (2 n\ 

^penguin 1 ^penguin 5 V / 



3 
N 

3,3.33 3 



A = diag(l ,1 ,1 , ), (2.4) 

Ajj $is$jd i 

where the first 3 (valence) entries of A are equal to 1 — 3/iV, and the next N + 3 (sea 
and ghost) entries are equal to —3/N. Here ip collects all quark fields in the theory, 
■0 = {QviiQshQgi)- The tensor A projects onto the valence (sd)i term in the first factor of 

Qpenguin- The motivation for Splitting Qpenguin this Way is that Qpenguin an d Qpenguin f° rm 

different representations of the partially quenched symmetry group: Q^nguin (Qpenguin) 
transforms in the trivial (adjoint) irrep of SU(3 + N\3)r. As a consequence, there are 
at least two different ways of embedding the QCD penguin operator into the partially 
quenched theory. One is to choose the partially quenched penguin to be a singlet under 
SU(3 + N\3)r, i.e. Qp e n gu iw as m the unquenched theory, whereas the other choice is to use 
the original operator, which is seen to be a linear combination of two irreducible operators. 
This latter choice was made in refs. [1, 2]. In ref. [10] non-singlet penguin operators such as 
Qpenguin were not considered, because singlet factors such as (uu+dd+ss)R in eq. (2.1) had 
been implicitly extended to singlets under the full (partially-) quenched symmetry group. 
Therefore, the analysis of ref. [10] was not complete, and ref. [3] and this paper remedy 
this for LR penguins. 

To leading order (LO), the operators Q^nguim Qp2i^uin are represented in ChPT by [3] 

QpeTguin - A) Stl (AL^) + Str (AX + ) , (2.5) 

QpelUn - f « (8 ' 8) str (AEAtf) , (2.6) 

where 

L u = iW^ , X± = 2B Q {Y I M^ ± MS f ) , (2.7) 
1 In a theory with K valence quarks, all ratios 3/N get replaced by K/N. 



3 



with M the quark-mass matrix, B the parameter B of ref. [11], E = exp(2«$//) the 
unitary field describing the partially quenched Goldstone-meson multiplet, and / the bare 
pion-decay constant normalized such that f n = 132 MeV. The a's are the correspond- 
ing LECs. Notice that Qp2^ U im unlike Qp2^ gu im * s °f order p°, due to the fact that the 
right-handed current in Qp2iguin * s n °t a partially quenched singlet (cf. electro- magnetic 
penguins 2 ). As already observed in ref. [3], the new operator Qp^ uin does not contribute 
at tree level to matrix elements with only valence quarks on external lines, since the matrix 
A is effectively proportional to the unit matrix in the valence sector. Indeed, replacing A 
by the unit matrix in eq. (2.6) would make the operator vanish. This is no longer true at 
next-to-leading-order (NLO), i.e. at order p 2 , where one-loop contributions from Qp2^ U in 
to valence-quark matrix elements are non-zero. 

Since the singlet-operator contributions also start at order p 2 , the NLO contributions from 
Qpenguin compete with the LO contributions from Qp2^ gu iw an d thus need to be taken into 
account already in a leading-order analysis of K — > 0, K — > it and K — > nn matrix elements. 
This also implies that a renormalization scale dependence already appears at leading order 
for those partially quenched matrix elements. That scale dependence is absorbed by new 
0(p 2 ) counterterms for Qp2^ U in- The complete list of C PS-even [4] operators is 

Qi QA = |^str(A{SASt,L,L M }) , (2.8) 



o(8,8) 

Q P 2 QA = str (AL^A^L,) , 



Qs QA = |^str(A{£,4£t,X + }), 



Ql QA = |— str(A[£A£t,X_]) , 



Ql QA = ^str(AEAEt) str (L M L /1 ) , 
Qe QA = |^str(ASASt)str(X + ) , 
Qr QA = ^i^str(A[EAEt,L M ]) , 
where we introduced the 0(p 2 ) LECs (3[ 8 ' 8 \. 

The partially quenched theory with N = 3 light sea quarks represents a special case. For 
N = 3, the LECs of the partially quenched theory must be the same as those of the 
physical, unquenched theory [12], basically because they represent the coefficients in an 



2 In fact, Qp2iguin i s a component of the same irrep as the electro-magnetic penguin, except for N = 

[3]- 
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expansion in powers of quark masses, and thus only depend on the number of dynamical 
(sea) quarks, and not on the quark masses themselves. Therefore, in the N = 3 partially 
quenched theory, what one should do is to omit Qp2^ U in altogether, because the aim is to 
obtain the values of af£\ and not the (8,8) LECs and /?f' 8) . 3 
The quenched case, N = 0, is different. In this case, the decomposition reads 



Q Q V en 9 u m = ^ str (A^ lu P L ) str (iP^Pr) + str (A^ 7 „P L ) str {N^ llx P R ) , 



1 

2 

= I nQ s + nQ NS (o q) 

2 ^penguin ' V penguin i V / 

N = ^diag(l, 1,1, -1,-1,-1) , (2.10) 

where valence entries of N are equal to |, and ghost entries are equal to — \. The first 
operator in the decomposition is a singlet under SU(3\3)r, while the second is not (NS 
for non-singlet). However, Qp^g Uin can transform into the singlet operator, implying that 
the non-singlet operators do not form a representation by themselves. In other words, 
Q penguin can m i x Ofpenguim which is possible because N is not supertrace-less, unlike A 
in the partially quenched case. The ChPT realization of both operators is obtained from 
the expressions given in eqs. (2.6,2.8) and by replacing A — > N. When referring to the 
quenched theory, we will add a subscript q to the LECs, and rename 4 the LECs for Qp^g Uin 

as a™ - <* A (8 ' 8) - P$ S - 

Within the fully quenched approximation, as also in the partially quenched case with 
N ^ 3, there is no reason that the LECs should have the same values as those of the 
unquenched theory. In general, the non-analytic terms are modified by quenching, and 
even the scale dependence of LECs is different between the quenched and unquenched 
theories. It is thus not a priori clear what choice to make for the embedding of penguin 
operators into the quenched theory. We will return to this issue in section 4 below. 



3 Non-singlet K° — > tt + tt matrix elements with general kine- 
matics 

In this section we present the partially quenched and quenched results for the contribution 
of strong penguin operators to the K° — > 7r + 7r~ matrix element to order p 2 in ChPT. We 
will restrict ourselves to the isospin limit in the valence sector, m vu = m v d, not assuming 
momentum conservation, so that q ^ Pi + P2, with q the (ingoing) K° momentum and p\ 
(p 2 ) the (outgoing) n + (7r~) momentum. All momenta are onshell and we work in euclidean 
space, i.e. p\= p\ = —M%, q 2 = —M\. M jvi (M jsi ) is the mass of a meson made out of 
the j'th valence quark and the ith valence (sea) quark. In all results presented below, we 
have symmetrized the expressions in the pion momenta p\ and pi- 

3 As long as one does not consider electro-magnetic penguin contributions. 

4 In the quenched theory there is no relation between Q£,^f u j„ and the electro-magnetic penguins [3]. 
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For the contribution of the singlet operator in eq. (2.5) a simple tree-level calculation yields 
the order p 2 result 

(n+n-\Q p Q s \K°) = —{a^ (±<z(pi + Pa) + P1P2) (3.1) 

M 2^ A , 1 2 Pl p 2 + q( Pl +p 2 )-2M 2 \ \ 
J V + 2 (g- Pl -p 2 ) 2 + Ml 

Note that the contribution proportional to a 2 vanishes for energy-conserving kinematics, 
because the corresponding operator can be written as a total derivative by using equations 
of motion [4] . As mentioned before, the tree-level contribution from the new non-singlet op- 
erator in eq. (2.6) vanishes. At one loop, it contributes at order p 2 . For general kinematics, 
we obtain, in the isospin limit, 

(n+n-lQ^lK )^ = ^[A(q;p u p 2 ) + A(q;p 2 , Pl )] , (3.2) 



A(q; Pl , P2 ) = ^a^Y,U(PiP2-M 2 )I(M 2 sl ,M 2 sl ,(p l+P2 f 

J i sea (. 

™l-Ml + m {L{ML) _ L(MLi 



(q-pif 
Ml 

(q-pi) 2 



-6 — + I f^ 9Pl + M ^ I(M 2 si , M 2 si , (q - Pl ) 2 ) 



(2 PlP2 - 2q Pl + Aqp 2 - 2M 2 ) 1 
+ (q - Pl - P2) 2 + M 2 K (L(MsJ L(M2J) | 
- (M jsi ^M jvi , j = 2, 3), 

i ghost 

where we used that M 2 si — M\ si = Ml — M 2 at this order, and 

/ 2 2 2 \ C ^ 

I(M u M 2 ,p) = - J ^ - M ^ _ -p - M ^ , 

where .D is the number of space-time dimensions. Explicit expressions for these integrals 
in the MS scheme are given in the appendix. The one-loop diagrams contributing to this 
result are displayed in fig. (1). In particular, the vacuum-tadpole diagram (last diagram 
in fig. (1)) gives a non-zero contribution for general kinematics which corresponds to the 
next-to-last line of eq. (3.2). In order to obtain the fully quenched result with no sea quarks 
at all, one simply drops the sum over sea quarks in eq. (3.2), keeping only the sum over 
ghost quarks, and replaces a^ 8 ' 8 ) — * a^ s . There are no contributions from rf double-poles. 

The new 0(p 2 ) counterterms in eq. (2.8) give, again for general kinematics 
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Figure 1: One-loop diagrams with the insertion of the non-singlet operator at the weak vertex 
(box) contributing to K° — > it + it~ with general kinematics 

_ V 8 > 8 )fM 2 _ M 2 ) + 2piP2 + q(pi+P2)-2M 2 \ 1 
3^ 3 ( ^ MJ i 2 {q-n-p.Y + Ml J/- (3 - 4) 

The structure of these contributions is such that the scale dependence contained in eq. (3.2) 
can be fully compensated by the LECs (3\ ' , i — 1, 2, 3. It is easy to verify that the other 
LECs (3\ 8,8 \ i = 4, . . . , 7 do not contribute. Since external legs only contain valence quarks, 
the tree-level K — > 7r7r matrix element does not involve the diagonal Ajj elements of the 
tensor A with % referring to a sea or ghost quark. This implies that, for this calculation, we 
may replace A with the unit matrix, and the operators Q p ® in eq. (2.5) for i — A, ... ,7 
vanish, while they become proportional to the operators in eq. (2.5) for % — 1,2, 3. It follows 
that the LECs af*^ together with fi[ 8 ^l will always appear in the specific combinations 

4 8,1) -(l-|)7^(2# 8) +^), (3.5) 

Q2 + V n) (4tt)^ 3 

for all tree-level matrix elements with only valence quarks on the external legs. For the 
analogue of eq. (3.5) in the quenched case, one replaces the factor 1 — 3/iV by 1/2, a^ 8 ' 1 ^ — > 
- < s and $ 8 ' 8) - 

We conclude this section with the expressions for the same matrix elements in the case of 
"physical" kinematics, i.e. with the choice q—pi + pi- Setting q = pi + P2 and using that 
q 2 = -M\ and p\=p\ = —M%, one obtains from eqs. (3.1), (3.2) and (3.4) 

(n + n-\Q p Q s \K°y h y s = j^af A) (M 2 K - M 2 ) , (3.6) 
(n + n-\Q PQA \KyX op = ^^^{4^ K I(M 2 si ,M 2 si ,-M 2 K ) (3.7) 

6 J i sea I Z 

^^-S^ I(M 2 Si ,M 2 Si ,-M 2 ) 



2 Ml 
'3 Ml 



2 M 2 



L(M 2 si ) - L(M 3 2 J 



- ]T (Mj s i >■ Mj v i, j = 2, 3) , 

i ghost 

{^-\Q p ^\ K y^ = (l-l) (2M 8 ' 8) +^ 8) ) (Ml - M 2 ) . (3.8) 
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Dropping the sea-quark terms and carrying out the sum over ghost quarks, we obtain a 
more explicit expression for the one-loop contribution in the quenched case, in the MS 
scheme: 

(n+n-lQ^lKX^ = (3.9) 

< / M 2 \ 



c*f 12(M|-M 2 )(log^-l 



16tt 2 / 3 q { \ A 2 

/Mt Mi ,\, M 2 

/Mt Mi 9 ,\, 2M 2 -M 2 

\ 7T 7T / 7T 

+2M^ (f(M 2 , M 2 , -Mi) - 2* vr 0(M| - 4M 2 ) ^1 - ^| + ^ 
+ - 2M|) (2F(M 2 , Mi, — M 2 ) + F(M|, 2M| - M 2 , -M 2 ))] 



where the function F is given in the appendix. 



4 Strategies for quenched estimates of real-world penguin ma- 
trix elements 

Recent numerical estimates of K — > nn matrix elements reported in refs. [1, 2] have been 
obtained via the indirect method, where the simpler K — > n and K — > transition ampli- 
tudes are computed on the lattice and then converted into estimates for K — > 7r7r matrix 
elements using ChPT. However, the fact that those numerical results are still obtained 
in the quenched approximation introduces a source of systematic error which is in princi- 
ple uncontrolled. As already explained, the LECs of the quenched theory do not have to 
have values equal to those of the unquenched theory. Typically, even the scale dependence 
of quenched and unquenched LECs is not the same; it depends on the number of light 
dynamical (sea) quarks in the theory. 

In the case of penguin operators, an additional ambiguity arises because it is a priori unclear 
whether it would be best to take matrix elements of Qp en f U i n , i-e. a linear combination of 
Qpenguin an d Qpmguin as m ec L- (2-9), or to drop the contribution from Q QNS under the 
assumption that a^ 1 ^ is the best estimate of a'j 8 ' 1 '*. 

In order to discuss possible strategies in more detail, we first recall the leading order ChPT 
expressions for K — > n and K — > matrix elements of strong penguin operators, from 
ref. [3]. In the quenched approximation (N = 0) one has 

{* + \qzuk + ) = AJ f - - j±y 2 + \f% s + ^) } , (4.1) 
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(0\Q? e C n ^n\K°) = j { (a^ + (Ml - Ml) (4.2) 

+< 5 E (HMD-HMD)} , 

i valence ) 

where Mk = M n = M in the case of the K — > 7r matrix element, and contributions of both 
singlet and non-singlet operators are included. Notice that only the (quenched versions 
of the) combinations (3.5) of LECs appear in these expressions, as expected. Assuming 
that one can limit the analysis to leading-order in ChPT, there are at least three different 
strategies for estimating K — > Tin penguin matrix elements from LECs obtained by fitting 
eqs. (4.1,4.2) to quenched numerical results: 

1. Ignore but not the other LECs associated with the non-singlet operator, /5^ 3 . 
Both /^g and (3 ql s + \(3q 2 S are sca le dependent (however, their sum is not, as can be 
seen from eq. (4.1)), implying that this strategy is scale dependent. However, it still 
makes sense in case the non-analytic contribution proportional to a^ s (cf. eq. (3.7)) 
is numerically small compared to all other contributions at a reasonable scale A of 
order 1 GeV. Thus, the linear combination — j^yi{f3qi S + |/9^ s ) is taken as the 

best estimate for the unquenched a\ , and eq. (3.6) can then be used to obtain the 
physical K — > -kit matrix element (at tree level). This is the strategy followed in 
refs. [1, 2]. In fact, in these works, it was assumed that the contribution proportional 
to a^ s in eq. (4.2) is small. 

2. Drop all the non-singlet operators. It was shown in ref. [3] that this can be done 
by dropping, in the fully quenched case, all eye-diagrams in which the right-handed 
quarks in eq. (2.1) are contracted. This can be easily deduced from eq. (2.9). This 
strategy was explored for Qq in ref. [13]. (For the partially-quenched case, see below.) 

3. Perform a complete quenched calculation including all contributions from singlet and 
non-singlet operators. After extracting all the LECs, singlet and non-singlet, one can 
use the sum of eqs. (3.6), (3.9) and the quenched version of eq. (3.8) to determine 
the quenched K — > -nit matrix element at the physical point. 

Strategy 2 isolates a q { , and might thus appear to be the obvious choice, since it is this 
LEC that is needed for calculating K — > tttt matrix element (to chiral leading-order) in 
the unquenched theory. However, as we already mentioned, the values of LECs in the 
quenched and unquenched theories do not have to be equal, and it might happen that 
(at some scale A) the quenched combination af^ — j^ysiftgi 3 + |/3^ s ), determined from 

strategy 1, is indeed a better estimate of af' 1 ^. Strategy 2 can be viewed as the situation 
in which the strong interactions are quenched at all scales between the weak and hadronic 
scales, because in that case only singlet penguin operators would appear in the evolution 
from the weak to the hadronic scale. So, while, on the one hand, it appears natural to 
assume only a mild flavor dependence of the LECs, in particular ai 8 ' 1 "*, one might, on the 
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other hand, argue that it is better to calculate the evolution from the weak to the hadronic 
scale in the unquenched theory, even if the matrix element at the hadronic scale is finally 
computed in the quenched approximation. The key point is that it is impossible to decide 
which strategy is best. 

The exception to these observations is the case of partially quenched QCD in which the 
number of light sea quarks is equal to that of the real world, in which N = 3. In the 
partially quenched theory, the singlet operator is (c/. eq. (2.3)) 

Qpe Q nguin = T7 Ml (u v ii v + d v d v + s v s v + ^ q si q si + u g u g + d g d g + s g s g \ , (4.3) 

V i ' R 

where the subscripts v, s and g denote valence, sea and ghost quarks, respectively. Strategy 
2 now corresponds to dropping all diagrams in which the right-handed valence and ghost 
quarks in the second factor of eq. (4.3) are contracted [3]. If the number of sea quarks 
N = 3 (but with the sea- and valence-quark masses not necessarily equal), the singlet 
LECs are those of the real world [12] , and therefore strategy 2 is the only correct one 
in this case. 

For any other case, fully quenched or partially quenched with N ^ 3, there is a priori no 
preferred choice; the spread in results obtained by employing all three strategies should 
be taken as a (lower bound of the) systematic error due to quenching. The extent to 
which strategies 1 and 3 lead to numerically different results depends on the size of 
contributions (at a given scale A). From eqs. (3.6,3.9) we find, taking physical values for 
all parameters, M K = 500 MeV, M n = 140 MeV, f = U = 132 MeV, M p = 770 MeV and 
M v = 550 MeV, that 

-i[K° -> 7T + 7T-], = 



400.7(c4i' 1) 


-af' 1 )) + (28.2- 


7.2i)< 5 


(A = 


1 GeV) 


400.7(c4 8 1 ' 1) 


-af' 1) ) + (32.2- 


7.2i)< 5 


(A = 


M P ) , 


400.7(c4 8 1 ' 1) 


-4 27 ' 1} ) + (37.3- 


7.2*)< 5 


(A = 





where we added in the tree-level ChPT contribution from the SU(3)l 27-plet operator [4]. 
If a^ s is of the same order as af{^ — af 1,l \ the contribution of the terms proportional to 
a™ is indeed small. The smallness of the coefficient of a^ s is due to a l/(47r) 2 suppression 
factor coming from the loop integral, and one might argue that /(47r) 2 is the "natural" 

parameter to compare with a^ 1 ^ — a^ 1,l \ in which case the contribution would not be 
small. Notice also that a small spurious imaginary part is generated by the non-singlet 
operator via the ghost-pion one-loop rescattering diagram. It is clear that the value of a^ s 
will have to be determined from a lattice computation. While this can be done by including 
the Oiq S terms, of e.g. eq. (4.2), in a fit to lattice data, there exists a much simpler and 
more reliable way of estimating the size of a^ s , as will be explained in the next section. 

Under the assumption that a^ s can be neglected without introducing a large uncertainty 
into the final estimate of strong penguin K — > tttt matrix elements, the question remains 
whether (to leading order in ChPT) or — -p^p- (Aji 5 + \ fiq2 S ) would be a better 
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estimate of af' 1 ^. The issue was investigated in ref. [13], where it was found that the 
difference between the two choices is numerically significant. At the physical kaon mass, 
the numerical value of the B parameter corresponding to Qq turns out to be approximately 
twice as large when the contribution of the non-singlet operator Qp^g Uin is omitted alto- 

(8 1) 

gether. Translated into estimates for the leading-order LECs, this implies that a q { is 
approximately twice as large as afi^ — {Pq\ S + \fiq2 S )- This may lead to substantial 
modifications in quenched estimates of e' /e, as discussed in ref. [5]. 

We emphasize that the whole discussion here is based on leading-order ChPT, and that 
NLO contributions may still lead to a substantial correction. However, it is reasonable to 
believe that NLO effects will not invalidate the basic content of our observations. 

Finally, we give a few more numerical examples of the partially quenched case with N = 2, 
always keeping the valence quark masses at their physical values (in the isospin limit), and 
choosing the two sea quarks to be degenerate in mass. Taking m sea = m u = rrid, we find 

-i[K° -> n + n-] N=2 = 400.7(af' 1) - a™) + 12.7«( 8 ' 8 ) (A = 1 GeV) , 

= 400.7(aS 8>1) - a (27 ' 1} ) + 14.0c/ 8 ' 8 ) (A = M p ) , 
= 400.7(a[ 8,1) - a™) + 15.7c/ 8 ' 8 ) (A = M„) , (4.5) 

whereas taking m sea = m s we obtain 

-i[K° -> 7r + 7r-]^ =2 = 400.7(af' 1) - a (27 ^) + (2.8 - 7.2i)a< 8 ' 8 > (A = 1 GeV) , 

= 400.7(ai 8,1) - a™) + (4.1 - 7.2i)a< 8 ' 8 > (A = M p ) , 
= 400.7(a£ 8,1) - a^ 7 ' 1 )) + (5.8 - 7.2*) a ( 8 ' 8 ) (A = M v ) . (4.6) 

Recall that for the N = 2 theory values of the LECs do not have to equal those of the N = 3 
theory. However, the partially quenched theory with two light sea quarks is closer to the 
real-world theory than the quenched (N = 0) theory. This is reflected by the fact that the 
coefficients of a*- 8,8 -* are small compared to those of a^ s in eq. (4.4). Notice in addition that 
in the N = 2 case, with m sea = m u = m d , the small spurious imaginary part vanishes, since 
it comes entirely from the pion-rescattering loop diagram where the sea-quark contribution 
is now fully cancelled by the corresponding ghost-quark contribution. 

In the case of N = 3 sea quarks with masses equal to the three valence quarks, ghost- and 
sea-quark contributions in eqs. (3.2,3.4) cancel, 5 ^[as they should, because this choice of 
parameters corresponds precisely to unquenched QCD. 



5 How to determine a™ on the lattice 

In principle, it is possible to determine a^ s from matrix elements with only physical 
(valence) particles as external states. For instance, given good enough statistics and a wide 

5 In cq. (3.4) the cancellation already occurs by just setting N = 3, because in this tree-level expression 
the sea- and ghost-quark masses do not appear. 
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enough range of quark masses, it can be determined from a fit to eq. (4.2). However, as also 
pointed out in ref. [2], the logarithmic terms in eq. (4.2) can look very linear in the typical 
range of quark masses used in lattice computations, making it hard to disentangle a^ s 

from af^ + (4^)2/^3 ^ ■ It would therefore be preferable to determine a^ s from a matrix 
element to which it contributes at order p°, because no other operators can "contaminate" 
the result at that order. 

It is very simple to do so, by considering matrix elements with ghost quarks on the external 
lines instead of valence quarks. Since this corresponds to a flavor rotation on the external 
lines, one needs to rotate the operator Qp^g Uin accordingly. A key point is that, while of 
course ghost quarks are not explicitly present in a quenched computation, their propagators 
are identical to those of the valence quarks, which are available in the actual computation. 

So, in order to determine a^ s , we propose to consider the following matrix element. First, 
we rotate Qp^g U i n by an SU(3\3)l rotation into 

QpeZguin = (s^PlS) Str (N^^Pr) . (5.1) 

This operator is in the same irrep of the group SU{3\3)l x SU(3\3) r , is thus parametrized 
by the same LECs as Qp^g Uin , and in particular, to leading order, by a^ s '. We then 

consider the matrix element of this operator between a fermionic kaon K oc ^755 and the 
vacuum. To leading order, 

mZUn\K) = 2ifa N q s + 0(p 2 ) , (5.2) 
thus isolating a^ s . Carrying out all quark Wick contractions, one finds that 

(0\Qinguin(y)d(xh 5 s(x)\0) = (5.3) 



l5(s(x)s(y)) lfli P L (d(y)d(x))\tT [ llx P R {{u{y)u{y)) + (d(y)d(y)) + (s{y)s(y)) 

+ {u{y)Ti{y)) + (d(yjd(y)) + (%)%)» 
-tr [^ 5 (s(x)s(y))-f^P R (s(y)s(y)}^^P L (d(y)d(x)} 
+tr [ 75 (s(x)s(y)) ^P L (d(y)d(y))^P R (d(y)d(x) ) 

where the traces are over spin and color indices only. A key observation is now that ghost 
propagators and valence propagators are equal flavor by flavor, (d(y)d(x)) = (d(y)d(x)) , 
etc.. Using this property, eq. (5.3) simplifies to 

(0\Q^ uin (y)d(xh 5 s(x)\0) = (5.4) 



-tr 



l5 (s(x)s(y)} lfl P L (d(y)d(x)}\ tr [7 ^P R ({u(y)u(y)) + (d(y)d(y)} + (s(y)s(y)}) 

1 



+ 2 tr 
1 

-2 tr 



-f 5 (s(x)s(y))^P R (s(y)s(y))^P L (d(y)d(x)) 
lh{s{x)s{y))^^P L (d{y)d{y))^^P R (d{y)d{x)) 



12 



We conclude that it is possible to estimate a^ s as a leading-order effect using only com 
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binations of contractions of valence-quark propagators. For the K — > matrix element 
of Qp^guin, cf. eqs. (2.9,4.2), the contractions in terms of valence quarks are of the same 
form, but the first two terms have the opposite sign, while the last term has the same sign 
as in eq. (5.4). Since the K — > matrix element is of order p 2 , we may combine the two 
results to obtain 



-tr 



75( S (^)%))7^L(rf(y)rf(y))7^(rf(l/)^))J ampMiated (5-5) 



= VZ ((OIQ^JK) + (0\Q^ uin \K)) 
= Vz(2 l fa» s + 0(p 2 )) , 



making it even easier to determine a^ s . The wave-function renormalization Z is defined 
from g?7 5 s = \J~Z~ K. 

The analysis for a similar determination of a^ 8 ' 8 ** in the partially quenched theory is anal- 
ogous. There, of course, the observation is not new, since a^ 8 ' 8 ^ is also the leading LEC 
for the electro-magnetic penguin (which in the partially quenched theory with N > 1 is in 
the same irrep as Qp2^ U in °f ec f (2-6) [3]). The main differences between a determination 
of o^ 8,8 ) and a^ s are that, first, a™ is not related to the electro-magnetic penguin in the 
quenched case [3], and second, that in order to determine it using leading-order (in this 
case 0(p°)) ChPT, one is forced to consider ghost quarks, as we did above. 



6 Conclusion 

In this paper, we continued our investigation of the ambiguities afflicting strong penguin 
contributions to K — > irn weak matrix elements due to the use of the quenched approxi- 
mation. 

The fact that the way of embedding penguin operators of the effective weak hamiltonian 
in the quenched theory is not unique tells us that, in the enlarged context of electro- weak 
interactions, the usual definition of the quenched theory is not complete. If only strong 
interactions are considered, it is sufficient to define quenched QCD as the modified version 
of QCD in which the quark determinant is set equal to a constant. A field-theoretic 
definition can be given through the introduction of ghost quarks into the path-integral 
[8], giving access to a complete picture of the symmetries of the quenched theory [9]. 
As soon as one considers operators external to QCD (i.e. the addition of electro-weak 
interactions), one has to answer the question how these operators should be incorporated 
into the quenched theory. Usually, this is straightforward. One classifies the operator 
by its flavor quantum numbers, in other words, one determines the irrep of SU(3)l x 
SU (3)a under which this operator transforms. If there exists a larger irrep of the quenched 
symmetry group which reduces to the unquenched irrep, the corresponding component of 
the quenched irrep can be taken as the quenched definition of the operator. However, 
in the case of strong penguins the operator, while irreducible in the unquenched theory, 
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is a linear combination of components of more than one irrep of the quenched symmetry 
group. Therefore each LEC of the unquenched theory corresponds to a set of LECs in 
the quenched theory. The ambiguity arises, because there is a priori no criterium for 
which linear combination of quenched LECs (if any) would yield the best estimate of 
the unquenched LEC. In the case of LR penguins considered here and in ref. [3], this 
phenomenon produces an effect already at leading order in ChPT. 

We remark that even in the simplest case when there exists a one-to-one correspondence 
between unquenched and quenched irreps, there is still the freedom to choose any compo- 
nent of the quenched irrep, and this flexibility can be used to extract LECs in the most 
convenient way [14]. However, in this case there is no ambiguity in the relation between 
unquenched and quenched LECs (even though their values may differ). This is in principle 
not different from the situation within the unquenched theory, where in general any com- 
ponent of an irrep can be used to extract the corresponding LEC. A classic example for 
weak matrix elements is the relation between B K and the K + — > 7r + 7r° decay rate [15]. (At 
non-leading order, it may not be possible to determine all LECs describing an operator in 
ChPT from one process, of course.) 

The ambiguity affecting penguin operators is fundamental, since there exists no solid theo- 
retical argument that can be used to decide the issue. Therefore, we argue that one should 
compare all choices that can be reasonably made, and take the resulting spread of estimated 
values as a lower bound on the systematic error due to quenching. It appears that in the 
case of e'/e this systematic error is rather large [13, 5]. A leading-order analysis of currently 
available lattice data [1, 2, 13] seems to indicate that quenched lattice computations cannot 
even confirm that this parameter is non-vanishing in the Standard Model. It could also 
be that the large numerical difference found between af^ and af^ — j^yz{(3qi S + \(3q\ S ) 
would be explained by the fact that higher orders in ChPT have not been taken into 
account, but we consider this to be unlikely. While it is clear that higher orders are nu- 
merically important, there appears to be no reason to assume that f3q\ s + \f3q 2 S is small. 
It could also be that a^ 5 , which appears in K° — > 0, and needs to be subtracted to obtain 

^i' 1 ^ — (4^)2 (@q\ s + \l3q\ s ) from K + — > 7r + , is not small. This would affect the determina- 
tion of a^ 1 ^ + j^yz(3q\ s , and hence the size of the subtraction. It is therefore important to 
obtain a reliable estimate of a^ s . We suggested a simple method for extracting its value. 

The above argument does not imply that lattice computations of e'/e are doomed to fail. 
On the contrary, quenched estimates of e'/ e with a particular choice for the strong penguins 
demonstrate that this computation is feasible, thanks to major advances in both theory 
and computational power. However, what will be needed in order to eliminate systematic 
errors due to quenching is a partially quenched study with N = 3 light sea quarks. This is 
the only approximation to unquenched QCD which is reliable in that it can be extrapolated 
systematically to the real world [12]. Currently existing quenched results give us invaluable 
information on what is needed to promote them to the required N = 3 world. For partially 
quenched QCD with N ^ 3, the situation is essentially the same as for quenched QCD, 
modulo differences in detail. 
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Summarizing, we presented the quenched and partially quenched results for the non-singlet 
contribution to the K° — > 7r + 7r~ matrix element. This made it possible to discuss in detail 
various strategies one might follow to use quenched computations in order to estimate 
the real-world value of this amplitude. Since a^ s contributes to this matrix element, but 
not to the K + — > 7r + transition amplitude used in ref. [1, 2], this introduces an additional 
ambiguity already at leading chiral order. The importance of this ambiguity depends on the 
size of oiq S and we have proposed a simple recipe for its determination. Our expressions for 
K° — > 7r + 7r~ with the most general possible kinematics and the inclusion of the non-singlet 
contributions, are appropriate for the analysis of direct quenched computations of this 
matrix element at leading order in ChPT. Beyond leading order, new problems arise [16], 
which may invalidate current methods for the direct determination of K — > -kit amplitudes 
with AI = 1/2 in quenched and partially quenched QCD. 



Appendix 



In this appendix we collect explicit expressions for the basic loop integrals appearing in 
eq. (3.2) etc. Using dimensional regularization, we have 



L(M 2 ) = J 



d D £ 



M 2 



(2tt) d £ 2 + M 2 16tt 2 
where A is the running scale, e = 4 — D, and 

Rei"(M 2 , M|,p 2 ) = -Re J -^L 



V 7 — log 4tt 

e 



, M 2 s 
+ log^-l 



(2n) D (i 2 + M 2 )((£ - p) 2 + M 2 ) 



(A.l) 



(A.2) 



16tt 2 \ 



+ 7 — log An 



M 2 1 / M\ M 2 \ M 2 1 



-1 + log-f + - 1 - + -f log-^ + -F(Mf,M 2 ,p 



A 2 



in which 

F(M 1 2 ,M|,p 2 



A 1, 



M[ Mf 

p2 ' p2 



log 



p 2 + M 2 + M 2 + p 2 ^/\(l,Mf/p 2 ,Ml/p 2 



p 2 + Ml + Ml - p 2 ^J\(l,Ml/p 2 ,Ml/p 2 ) ' 
X(x,y,z) = (x - y + zf + Axy . (A.3) 



MS expressions are obtained by dropping the contact terms in square brackets. 

For p 2 > 0, the argument of the logarithm in eq. (A.3) is positive, and I is real. For 
P 2 < 0, F is obtained by analytic continuation. A(l, Ml/p 2 , M^/p 2 ) turns negative for 
-(Mi + M 2 ) 2 < p 2 < -(Mi - M 2 ) 2 , and we find that I is still real with F now given by 



F(Mf,M 2 ,p 2 ) = 2 



. / Ml M 2 \ -p 2 ^/-\(l,Ml/p 2 ,M 2 /p 2 ) 
" A I h 1? ) arCtan j? + M? + M$ 



(A.4) 
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At p 2 = —{Ml + Mf) the argument of the arctangent has a singularity, across which the 
branch of the arctangent has to be chosen continuously: 



arctan • 



-p- 



-\(l,M?/p*,M%/p>) 



p 2 + Ml + Ml 



Arctan 



-p 2 ^/-\(l,Ml/p 2 ,M 2 /p 2 
p 2 + Ml + Mf 
P 2 < -{Ml + M 2 2 ) , 



+ 7T , 



(A.5) 



where Arctan denotes the principal value of the arctangent. Again continuing analytically 
across p 2 = — {M 1 + M 2 ) 2 , F is again given by eq. (A. 3), but 7(M 2 ,Mf,p 2 ) picks up an 
imaginary part: 



Im I {Ml, M 2 ,p 2 ) = 



16n 2n \ 



A 1 



Ml Ml 



(A.6) 
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